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Optimal  Allocation  of  Multistate  Components 

by 

Emad  El-Neweihi,  Frank  Proschan  and  Jayaram  Sethuraman 


ABSTRACT 


In  this  paper  we  present  some  results  in  the  optimal  allocation  of 
multistate  components  to  k  series  systems  so  that  some  performance  charac¬ 
teristic  like  expected  number  of  systems  functioning  at  level  a  or  higher, 
the  probability  that  at  least  one  of  the  systems  functions  at  level  a  or 
higher,  etc.  is  maximized.  Our  basic  mathematical  tools  are  majorization 
and  Schur  functions;  the  methods  used  and  some  of  the  theorems  obtained  are 
those  of  "Optimal  Allocation  of  Components  in  Parallel-Series  and  Series- 
Parallel  Systems"  El-Neweihi,  E.  ,  Proschan,  F.,  and  Sethuraman,  J.,  Report 
(1984).  In  addition,  we  show  how  these  results  may  be  used  to  obtain  fruit¬ 
ful  applications  in  reliability  theory. 


1.  Introduction. 


The  theory  of  multistate  systems  has  been  developed  to  describe  more 
adequately  the  performance  of  systems  and  their  components  each  of  which 
operates  at  more  than  the  usually  assumed  two  levels  of  performance: 
"functioning"  or  "failed".  (For  a  recent  survey  of  the  subject  see  [3].) 

In  this  paper  we  present  some  results  in  the  optimal  allocation  of 
multistate  components  to  k  series  systems  so  that  some  performance  charac¬ 
teristic  like  expected  number  of  systems  functioning  at  level  a  or  higher, 
the  probability  that  at  least  one  of  the  systems  functions  at  level  a  or 
higher,  etc.,  is  maximized.  Our  basic  mathematical  tools  are  majorization 
and  Schur  functions;  the  methods  used  and  some  of  the  theorems  obtained  are 
those  of  [2].  In  addition,  we  show  how  these  results  may  be  used  to  obtain 
useful  applications  in  reliability  theory.  The  present  paper  and  reference 
[2]  are  the  only  papers  exploiting  the  elegance  and  power  of  majorization  and 
Schur  functions  to  solve  optimal  allocation  problems  in  reliability  as  far 
as  we  know. 

Preliminaries . 

For  vector  x  =  (x^,  ...,  xn)  ,  let  x^-j  2  ...  2  Xj-^-j  denote  the  decreasing 
rearrangement  of  the  coordinates  of  x.  We  say  vector  x  majorizes  vector  y  if 

k  k 

[  x[i]  2  I  ’’[i]  f°r  k  *  ••••"-  ! 

and 

n  n 

\  X[il  =  £  y[i]' 

We  say  a  function  f:  Rn  -*•  R  is  Schur-convex  (Schur- concave)  if 
f(x)  ^  f(y)  (f(x)  £  f(y))  whenever  x  5  y , 


Examples  of  Schur-convex  (Schur- con cave)  functions  are  £  fj(x-)  and 
n  i  =  1  1  1 

XJ  f2(x^),  where  f ^ :  R  -*■  R  is  convex  (concave)  and  f ^  ■  R+  R  is  log-convex 
i  =  1 

(log- concave) . 

A  random  variable  X  is  said  to  be  stochastically  larger  than  a  random 

st 

variable  Y  if  P(X>x)  2  P(Y  >  x)  for  all  real  x;  we  write  X  >  Y.  For  vectors, 

X  ^  Y  if  f(X)  f(Y)  for  every  increasing  function  f. 

A  multistate  series  system  consisting  of  n  raultistate  components  is  a 

system  whose  structure  function  is  given  by  <j>(x)  =  rain  x.  where  x.  represents 

l<i<n  1  1 

the  state  of  component  i  and  takes  its  value  in  a  common  state  space  S  £  [0,») 
i=l,  ...,  n.  A  multistate  parallel  system  can  be  similarly  defined. 


2 .  Optimum  Allocation  of  Multistate  Components. 

We  review  a  general  optimal  allocation  result  for  multistate  systems 
recently  obtained  in  [2].  We  also  present  additional  related  results.  We 
then  describe  in  some  detail  reliability  models  in  which  these  results  can  be 
used. 

The  following  theorem  is  basic  in  [2];  the  proof  depends  on  the  useful 
tools  of  majorization  and  Schur  functions.  We  present  an  outline  of  this 
proof;  the  interested  reader  may  consult  [2]  for  further  details. 

2,1.  Theorem.  Let  P^ ,  _ ,  P^  be  the  disjoint  min  path  sets  of  a  parallel- 

series  system  having  path  lengths  n^ ,  ...,  n^.  Without  loss  of  generality, 
assume  that  n^  <  ...  <  n^.  Suppose  that  there  are  n  =  ^  ^  +  ...  +  n^ 

independent  components  with  reliabilities  p^ ,  ...,  p^  (at  time  tQ,  say)  to 
be  allocated  among  the  path  sets.  Then  the  reliability  of  the  system  (at 
time  t  )  is  maximized  when  the  n^  most  reliable  components  are  allocated 
to  P^,  the  n^  next  most  reliable  components  are  allocated  to  P-,,  ...,  and 
finally,  the  n,  least  reliable  components  are  allocated  to  P,  . 


Proof.  Let  x,.  =  log  JT  p.,  i  =  1,  ....  k.  The  allocation  described  in  the 
1  j  e  P.  J 

hypothesis  maximizes  x=  (x  ,  x,)  in  the  sense  of  majorization.  The  reli- 

i  k  k 

ability  function  of  the  system,  1  -  JY  (l-e-xi)  is  a  Schur-convex  function 

i=l 

of  x,  and  hence  the  reliability  is  maximized  by  the  allocation  described. || 

The  following  theorem  shows  that  the  allocation  described  above  also 
maximizes  the  expected  number  of  working  min  path  sets  which  are  now  viewed 
as  k  separate  series  systems. 

2.2.  Theorem.  Let  S^,  . ..,  be  k  series  systems  of  sizes  n^  <  ...  ^  n^, 
respectively.  Suppose  that  there  are  n  =  n^  +  ...  +  n^  independent  components 
with  reliabilities  p^,  . ..,  pn  (at  time  tQ,  say)  to  be  allocated  among  the 
series  systems.  Then  the  expected  number  of  working  subsystems  (at  time  tQ) 
is  maximized  when  the  n^  most  reliable  components  are  allocated  to  S^,  the 
n.,  next  most  reliable  components  are  allocated  to  S?,  . ..,  and  finally,  the 
n^  least  reliable  components  are  allocated  to  S^. 

Proof.  Let  x  =  (x, ,  . ..,  xu)  be  as  defined  in  the  proof  of  Theorem  2.1.  The 

“  k  x. 

r  i 

expected  number  of  working  systems  is  [  e  which  is  a  Schur-convex  function 

i  =  1 

of  )c.  The  result  thus  follows,  jj 

Note  that  the  optimal  allocation  in  Theorem  2.1  and  Theorem  2.2  does  not 
depend  on  the  actual  values  of  the  reliabilities  p^,  ...»  pn>  but  only  on 
their  ordering. 

The  above  results  can  be  extended  to  cover  the  case  of  multistate  components 
and  systems.  Consider  a  parallel-series  system  as  described  in  Theorem  2.1, 
except  that  now  each  component  has  a  common  state  space  Sc  [0,®);  common 
choices  for  S  are  (0,  1,  . . . ,  M)  and  the  unit  interval  [0,1].  Let  X^(t), 

...,  Xn(t)  denote  the  states  of  the  components  1,  ...,  n  and  X(t)  denote  the 
state  of  the  system  at  time  t,  t  >  0.  Suppose  that 
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(2.1)  X1(t)  <  X2(t)  <  ...  V  Xn(t) 

for  each  t  >  0.  Theorem  2.3  just  below  describes  the  best  allocation  of  com¬ 
ponents  to  path  sets  to  maximize  the  system  state  X(t)  in  the  stochastic  sense 
uniformly  in  t. 

2 , 3.  Theorem.  Consider  a  parallel-series  system  as  described  above  with  com¬ 
ponent  states  satisfying  (2.1).  To  maximize  the  system  state  stochastically, 
uniformly  in  t,  allocate  components  as  follows:  The  best  n^  components  (those 
with  stochastically  largest  states)  to  P  ,  the  next  best  n2  components  to 
P7,  .  ..,  and  finally  the  n^  worst  components  to  P^. 

Proof.  Fix  a  time  t  >  0  and  state  a  2:  0.  Let  P“(t)  =  PfX^t)  >  a) , 

1  <  i  <  n.  These  probabilities  are  in  increasing  order  in  i  for  any  fixed 
t  and  a  in  view  of  condition  (2.1).  Since  the  optimal  allocation  in  Theorem 
2.1  depends  only  on  the  ordering  among  such  probabilities,  the  same  allo¬ 
cation  maximizes  the  system  state  stochastically,  uniformly  in  t.  || 

Now  let  us  view  the  above  min  path  sets  as  k  separate  multistate  series 
systems  and  call  them  S^,  ...,  S^.  For  ct  2  0,  let  Na(t)  be  the  number  of 
systems  functioning  at  a  level  a  or  higher.  We  have  the  following  theorem. 


2.4.  Theorem.  Let  ...,  SR  and  Na(t)  be  as  described  above.  Then  the 

allocation  described  in  Theorem  2.3  maximizes  ENa(t)  uniformly  in  a  and  t. 

k 

Proof  Fix  a  time  t  >  0  and  a  state  a  2  0.  Then  E(Na(t))  =  £  IT  P-(t), 
-  i-1  jeSi  J 

where  P®(t)  =  P(Xj(t)  >  a)  ,  j  =  1 ,  • • • ,  n.  The  result  now  follows  from  Theorem 
2.2  and  Condition  (2.1).  || 

In  the  results  we  have  considered  so  far,  all  the  components  are  tacitly 
assumed  to  be  of  the  same  type  and  thus  can  be  interchanged  freely.  We  now 
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consider  a  situation  in  which  we  have  multistate  components  of  m  different 
types.  Let  xj(t)  s<^  X^(t)  ...  s<?  (t)  be  the  states  of  k  components  of 

type  j  at  time  t,  j  =  l,  i,  t  >  0.  We  wish  to  build  k  multistate  series 

systems  each  containing  one  component  of  each  type.  Let  Na(t)  be  the  number 
of  systems  that  function  at  level  ct  or  higher,  where  a  >  0  corresponds  to  an 
arbitrary  but  fixed  state.  We  wish  to  maximize  either  E(Na(t))  or  P(Na(t)  2  1). 
The  following  theorem  gives  the  optimal  allocation. 

2.5.  Theorem.  Let  S^,  ....  denote  the  k  series  systems.  To  maximize 
E(Na(t))  or  P(Na(t)  2  1)  uniformly  in  a  and  t,  assemble  the  systems  as 
follows : 

The  best  components  are  assembled  together,  the  next  best  components 
are  assembled  together,  and  finally  the  worst  components  are  assembled 

together.  Thus  system  is  built  with  components  whose  states  are  xf(t), 

....  X™(t],  i-1,  ....  k. 

Proof.  Let  P? ,a(t)  =  P(x|(t)  2  a)  ,  j  =  1,  . . . ,  m  and  i =  1,  . . . ,  k.  Let 

m  •  j 

xa  =  log  HV?01  ,  i  =  1 ,  . . . ,  k,  where  <p  ,  . . .  ,  are  m  permutations  of 

1  j  =  l  *J(i) 

{1,  ...,k}.  Let  x3  =  (x^,  x£).  The  assembly  described  in  the  hypothesis 

of  the  theorem  maximizes  xa  in  the  sense  of  majorization .  Since  both  E(Na(t)) 
and  P(Na(t)  2  1)  can  be  expressed  as  Schur-convex  functions  ofxa,  the  result 
now  follows.  || 

2.6.  Remark .  The  model  discussed  just  above  is  related  to  a  binary  model 
given  in  [1]  . 

3.  Applications  in  Reliability. 

In  this  section  we  describe  models  to  which  theorems  of  section  2  can 
be  applied.  In  all  models  the  state  space  is  {0,  1,  ....  M}. 
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3.1.  Models.  Consider  n  independent  binary  components.  Suppose  that  each 
component  is  supported  b>  M  -  1  functioning  spares  that  do  not  deteriorate 
until  put  into  use.  We  say  that  "position"  or  "socket"  i  is  in  state  M  when 
the  original  component  of  type  i  is  still  functioning  and  none  of  the  M-  1 
spares  have  been  used;  the  position  is  in  state  M-  1  when  the  original  com¬ 
ponent  has  failed  and  has  been  replaced  by  a  spare,  leaving  M-  2  spares 
available  for  replacement;  ...;  and  finally,  the  position  is  in  state  0  when 

no  spares  are  available  and  the  component  in  use  fails. 

Next  let  be  the  life  length  of  the  spare  for  component  type  i, 

1  <  i  <  n,  1  <  j  <  m.  Note  that  each  original  component  in  operation  at 
time  0  is  viewed  as  a  member  of  its  spares  kit.  Assume  that  for  each  i, 

T*  .  .  .  ,  T*  are  independently  and  exponentially  distributed  with  parameters 
A*,  ...,  A*  respectively.  Let  Xtl)  -  X(2)  -  •  -  X(M)  be  a  rearrangement  of 

a|,  . . .  ,  Ajjj  ,  i  =  1,  . . .  ,  n.  Assume  that  A^  >  A^  >  ...  >  >  ^  for  l-  1,  .  . .  ,  M. 

Note  that  for  each  pair  (i,t),  the  distribution  of  X^(t)  depends  on  the 
order  in  which  we  have  been  using  the  spares.  However,  a  little  reflection 
shows  that  the  following  order  maximizes  X^(t)  stochastically,  uniformly  in  t 
for  each  i:  Start  with  the  spare  whose  parameter  is  A^  .  Upon  its  failure, 
replace  it  with  the  spare  whose  parameter  is  A^,  ...,  and  finally  use  the 

spare  whose  parameter  is  A^.  Let  Xt(t)  be  the  random  state  of  component  i 

corresponding  to  the  above  order.  Then  clearly  X*(t)  S<"  X*(t)  S<"  ...  S<"  X*(t) 
for  each  t.  The  allocation  described  in  Theorem  2.3  maximizes  stochastically 
the  system  state  uniformly  in  t. 

Thus,  employing  the  order  described  above  for  using  the  spares  corres¬ 
ponding  to  each  socket,  together  with  the  optimal  allocation  of  Theorem  2.3 
yields  the  best  stochastic  performance  of  a  parallel-series  system  formed 
from  the  n  multistate  components.  [| 
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Similarly  Theorem  2.4  gives  us  the  allocation  that  maximizes  the  expected 
number  of  Na(t)  (uniformly  in  a  and  t) ,  where  Na(t)  is  the  number  of  systems 
functioning  at  a  level  a  or  higher  among  a  total  of  k  series  systems  built 
from  the  above  multistate  components. 

3.2.  Remark .  In  the  example  given  in  [2]  ,  X*  =  X*  =  . . .  =  X^  =  X^ ,  i  =  1 ,  . . .  ,  n. 
Thus,  in  this  case,  the  order  in  which  the  spares  are  used  is  immaterial. 

3.3.  Remark .  In  Model  3.1,  it  is  assumed  that  for  each  i,  T*,  . . .  ,  T*  have 
exponential  distributions  with  parameters  satisfying  certain  conditions.  Other 
distributions  can  be  used  provided  the  following  conditions  hold:  1)  For  each 
i,  the  random  variables  T^,  ...,  T^j  are  stochastically  ordered.  2)  For  i<  j, 
the  smallest  stochastically  among  T*,  ...,  T*  is  stochastically  less  than  the 
smallest  among  t|,  ...,  T^,  the  second  smallest  among  T*,  ...,  T*  is  stochas¬ 
tically  less  than  the  second  smallest  among  t|,  ...,  T^,  ....  and  finally  the 
largest  stochastically  among  T^,  . . . ,  Tjjj  is  stochastically  less  than  the 
largest  stochastically  among  T|,  ...,  T^  . 

3.4.  Model.  Consider  n  •  M  binary  components  forming  n  parallel  systems 

Sj,  . . . ,  Sn  of  size  M  each.  We  now  view  each  as  a  single  multistate  system 
of  components  with  M+  1  states  0,  1,  ...,  M,  as  follows:  When  all  the  M 
binary  components  in  system  are  functioning,  then  the  socket  corresponding 
to  multistate  component  i  is  in  state  M.  When  the  first  binary  component  in 
system  fails,  socket  i  is  now  in  state  M-  1,  and  so  on,  until  the  last 
binary  component  in  system  fails,  and  socket  i  is  now  in  state  0,  i=  1,  ...,  n. 
Let  T 1  =  (T*,  ...,  Tjjp  be  the  random  vector  representing  the  joint  lifelengths 
of  the  M  binary  components  in  system  S^,  i=  1,  ...,  n.  Suppose  that 
T1  S<"  T^  s<^  ...  S<"  Tn  .  Let  X^(t),  ...,  X  (t)  be  the  states  of  components 
1,  . . . ,  n  at  time  t,  t  >  0.  Then  P(X^(t)  >  j)  =  P(T^  ^  >  t)  for  each  j  =  0, 
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M  —  1 ;  i  =  1,  n,  and  t  ^  0  where  is  the  l- th  order  statistic 

among  T* ,  .  ..,  Tjjj.  Clearly  X^(t)  S<?  X2U)  S<?"  ...  s«?  X^ft)  for  each  t  ^  0. 

Theorem  2.3  can  now  be  used  to  provide  us  with  the  optimal  allocation  of  the 
n  multistate  components  to  the  min  path  sets  of  a  parallel-series  system  of 
n  components.  || 

3,5.  Model.  Consider  k  binary  components  of  type  j,  j  =  1,  m.  The 

i*-*1  component  of  type  j  has  a  spares  kit  with  M  -  1  functioning  spares  that 
do  not  deteriorate  until  put  into  use,  i=l,  ...,  k,  j  =  1,  . . . ,  m.  We  view 
each  binary  component  together  with  its  spares  as  a  multistate  component  in 
the  same  fashion  as  described  in  Model  3.1.  The  life  length  of  the  i^  component 
of  the  type  as  well  as  the  life  lengths  of  all  of  its  spares  have  exponen¬ 
tial  distributions  with  the  same  parameter  X?,  i=l,  ...,  k;  j  =  l,  ...,  m.  All 
the  random  variables  considered  are  independent.  Assume  x|  £  X^  >  . . . >  X^  , 
j  =  1,  . . .  ,  m.  Now  let  X^(t)  be  the  state  of  the  ith  component  of  type  i  at 
time  t,  where  i  =  1,  k,  j  =  l,  ...,  m,  t^O.  Clearly  x|(t)  S<"  X^(t)  s<" 

. .  .  Xj^(t) .  Theorem  2.5  can  now  be  used  in  the  obvious  fashion  to  obtain  the 
optimal  assembly  of  k  series  system.  || 
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